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Abstract. We consider the propagation of wave packets for the nonlinear 
Schrodinger equation, in the semi-classical limit. We establish the existence 
of a critical size for the initial data, in terms of the Planck constant: if the 
initial data are too small, the nonlinearity is negligible up to the Ehrenfest 
time. If the initial data have the critical size, then at leading order the wave 
function propagates like a coherent state whose envelope is given by a nonlinear 
equation, up to a time of the same order as the Ehrenfest time. We also prove 
a nonlinear superposition principle for these nonlinear wave packets. 



1. Introduction 

We consider semi-classical limit e — > for the nonlinear Schrodinger equation 
iied t r + Y A ^ =V{xW + \\r\ 2 °r, (M)eR+xR d , 

with A G R, d ^ 1. The nonlinearity is energy subcritical (a < 2/{d — 2) if d ^ 3). 
This equation arises for instance as a model for Bose-Einstein Condensation, where, 
among other possibilities, V may be exactly a harmonic potential, or a truncated 
harmonic potential (hence not exactly quadratic); see e.g. [13], I26j. 

Assuming that V is smooth and subquadratic (this notion is made precise below, 
see Assumption ll.ip . we know that for each e > 0, (jl.ip has a unique global solution 
in the energy space 

S={/eH 1 (R d ), x^\x\f(x)eL 2 (R d )} 1 

provided -0o G S and, either a < 2/d, or (er ^ 2/d and A ^ 0), while if A 2/d 
and A < 0, finite time blow-up may occur; see [7J. We assume that the initial data 
ipQ is a localized wave packet of the form 

(1.2) Vo(z) = £ P x £~ d/4 a (^p) e i ^- Xo >^ a , a £ S(R d ), x 0) £o G R d . 

Such data, which are called semi- classical wave packets (or coherent states), have 
been extensively studied in the linear case (see e.g. [4j [TU1 [Til [33j [35] ) . I n particular, 
Gaussian wave packets are used in numerical simulation of quantum chemistry like 
Initial Value Representations methods. On this subject, the reader can refer to 
the recent papers [3S1 [37J [3S] where overview and references on the topics can be 
found. These methods rely on the fact that if the data is a wave packet, then the 
solution of the linear equation (A = 0) associated with (jl.l[) still is a wave packet at 
leading order up to times of order Clog (=■): such a large (as e — > 0) time is called 
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Ehrenfest time, see e.g. [51 122} [23] ■ Our aim here is to investigate what remains of 
these facts in the nonlinear case (A ^ 0). 

In the present nonlinear setting, a new parameter has to be considered: the size 
of the initial data, hence the factor e@ in (|1.2p . The goal of this paper is to justify a 
notion of criticality for (3: for (3 > (3 C :— l/(2a) + d/4, the initial data are too small 
to ignite the nonlinearity at leading order, and the leading order behavior of ip 6 &s 
e — ► is the same as in the linear case A = 0, up to Ehrenfest time. On the other 
hand, if /? = /3 C , the function ip £ is given at leading order by a wave packet whose 
envelope satisfies a nonlinear equation, up to a nonlinear analogue of the Ehrenfest 
time. We show moreover a nonlinear superposition principle: when the initial data 
is the sum of two wave packets of the form (|1.2|) . then ip e is approximated at 
leading order by the sum of the approximations obtained in the case of a single 
initial coherent state. 

Up to changing ip £ to e~^i\) e , we may assume that the initial data are of order 
0(1) in L 2 (R d ), and we consider 



(1.3) 



F 2 

ied t ip E + — A^ £ = V{x)ip e + Ae Q |i/; e | 2<T V e , (t, x) G R+ x R d , 
^ e (0,x) =£- d / 4 a (^p) e^-^'W*, 



where a = 2/3(7 . 

1.1. The linear case. In this paragraph, we assume A = 0: 

£2 A./.e rr/.\./. E . „s _ p -d/i n ( X-Xq \ A( x - X0 ).£ Q / £ 



(1.4) ied t iP e + — A^ e = V(x)ip e ; V e (0, x) = £~ d/4 a J e 

The assumption we make on the external potential throughout this paper (even 
when A 0) is the following: 

Assumption 1.1. The external potential V is smooth, real-valued, and subquadratic: 
V G C°° (R d ; R) and d^V £ L°° (R d ) , V| 7 | ^ 2. 

If I 2 

Consider the classical trajectories associated with the Hamiltonian ^ — h V(a;): 

(1.5) x(t) = £(t), ffl = -W(i(t)); i(0) = io, C(0) - Co- 
These trajectories satisfy 

l -^ + V(x(t)) = l *f + V(x ), ViGR. 

The fact that the potential is subquadratic implies that the trajectories grow at 
most exponentially in time. 

Notation. For two positive numbers a e and 6 e , the notation a e < b e means that 
there exists C > independent of e such that for all e G]0, 1], a £ ^ Cb £ . 

Lemma 1.2. Lei (a^oj^o) G R d x R d . Under Assumvtion M . 1[ (|1.5|) /ias a unique 
global, smooth solution (x, £) € C^RjR^) 2 . 7i grows at most exponentially: 

(1.6) 3C >0, |ar(t)| + <c Cot , VieR. 
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Sketch of the proof. We explain the exponential control only. We infer from (|1.5|) 
that x solves an Hamiltonian ordinary differential equation, 

x{t) + W (x(t)) = 0. 

Multiply this equation by x(t), 

l((if + V(x(t)))=0, 

and notice that in view of Assumption ll.li V(y) < (y) 2 '- 

m < (x( t )) , 

and the estimate follows. □ 
Remark 1.3. The case V(x) = —\x\ 2 shows that the result of Lemma ITT21 is sharp. 
We associate with these trajectories the classical action 

(1.7) S(t)=J*(±\Z(s)\ 2 -V(x(s))) ds. 

We observe that if we change the unknown function tjj £ to u e by 

(1.8) ^(t, x) = £- d/ V (t, ^=^) e *(SW+CW-(x-x(t)))/= j 
then, in terms of u e — u E (t,y), (|1.4p is equivalent 

(1.9) l c» tU £ + iA u e = y £ (t,zjK ; " £ (0,y) =a(tf), 
where the external time-dependent potential y e is given by 

(1.10) V e (t, y) = - £ (V(x(t) + y/iy) - V{x{t)) - yfe (VV(x(t)), y)) . 

This expression reveals the first terms of the Taylor expansion of V about the 
point x(t). Passing formally to the limit, V e converges to the Hessian of V at x(t) 
evaluated at (y, y). One does not even need to pass to the limit if V is a polynomial 
of degree at most two: in that case, we see that the solution ip e remains exactly a 
coherent state for all time. Let us denote by Q(t) the symmetric matrix 

Q(t) = HessV(x(t)). 

It is well-known that if v solves 

(1.11) id t v + ^Av = ^{Q{t)y,y)v ; v(0,y) = a{y), 
then the function 

(1.12) ^f in (i, x) = e- d ' A v [t, ^|^) e «s<t)+«t).(«-x(t)))/ e 

approximates tp e for large time in the sense that there exists C > independent of 
e such that 

\\r(t,-)-V>Ut,-)\\L H n*)<CV^e Ct . 

See e.g. EOl HU HH HU [2H HH] and references therein. We give a short proof 
of this estimate, which can be considered as the initial step toward the nonlinear 
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analysis which is presented in the next paragraph. We first notice that since V is 
subquadratic, we have the following pointwise estimate: 



(1.13) 



V e (t,y)-~(Q(t)y,y) 



< QVZ\y\ 



for some constant C_ independent of t. The error rf in = u e — v satisfies 

id t rl n + \ Arf in = V s u s - \ (Q(t)y, y) v = V e r^ + (V - \ (Q(t)y, y)) v, 

along with the initial value rf in i t=0 = 0. Since V £ is real-valued, the classical energy 
estimate for Schrodinger equations yields, in view of (jl.131) . 

lkfmlU-([o,t];^(R^)) < Ve / W\y\ 3 v(T,y)\\ L 2 {Ild) dT. 

Since Q is bounded (V is subquadratic), we have the control 

\\\y\ 3 v(T,y)h*(n«)<Ce CT 

for some constant C > 0; see Proposition 12. II below. We then have to notice that 
the wave packet scaling is L 2 -unitary: 

writ, •) - vL(t, oii^cr*) = n« £ (^ ■) - «(*» oiil^r")- 

To summarize, we have: 

Lemma 1.4. Let d 1 cmd a G 5(R <i ). There exists C > independent of e such 
that 

(1.14) ||^(t, •) - pf in (i, Oll^ca^) < C^e c *. 
In particular, there exists c > independent of e such that 

sup ||V> £ (t, •) - tpf in (t, -)IU=(R<*) — >0. 

1.2. The nonlinear case. We now consider the nonlinear situation A ^ 0. Re- 
suming the same change of unknown function (|1.8|) , then adapting the above com- 
putation leads to 

(1.15) id t u £ + ^Au £ = V £ u e + \s a ~ a < \u e \ 2a u e , 

where V e is given by (|1.10[) as in the linear case, and 

da 

(1.16) a c = l + — . 

The real number a c appears as a critical exponent. In the case a > a c , we can 
approximate the nonlinear solution u e by the same function v as in the linear case, 
given by (jl.lip . The space E will turn out to be quite natural for energy estimates. 
Introduce the operators 

A £ (t) = V-eV-M ; B £ (t)= X -p±. 

Note that A and B are essentially V and x, up to the wave packet scaling, in the 
moving frame. From this point of view, our energy space is quite different from the 
one associated with the Lyapounov functional considered in [2] , and more related 
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to the one considered in |25j , since we pay attention to the localization of the wave 
packet, through B £ . For / 6 S, we set 

WfWn = + \\A £ (t)f\\ L2{Ild) + \\B £ (t)f\\ LHRd) , 

a notation where we do not emphasize the fact that this norm depends on e and t. 

Proposition 1.5. Let d ^ 1, a G S(R d ). Suppose that a > a c . There exist 
C, C\ > independent of e, and Eq > such that for all e £E]0, £o], 

\\ip £ (t) - ff in (t)\\ n < e 7 e Clt , O^ts^Clog-, where 7 = min [ i a - a, 

V- 

In particular, there exists c > independent of e such that 

sup \\r(t)-vL(t)\\n— >0- 

The proof is more complicated than in the linear case (see Q . The solution of 
(|1.3[) is linearizable in the sense of [TB] (see also [H]), up to an Ehrenfest time. 

In the critical case a = a c with A 7^ 0, the solution of (|1.3|) is no longer lineariz- 
able. Indeed, passing formally to the limit e — > 0, Equation (|1.15[) becomes 

(1.17) id t u+^Au = ^{Q(t)y,y)u + \\u\ 2 °u ; u(0, y) = a(y). 

Remark 1.6 (Complete integrability) . The cubic one-dimensional case d = a = 1 is 
special: if Q = 0, then (|1.17|) is completely integrable ([I]). However, if Q ^ 0, there 
exists no Lax pair when the nonlinearity is autonomous as in (|1.17|) ; see [3UJ [35J • 
Note also that if Q = 0, then u £ = u for all time. 

As in the linear case, we note that if V is exactly a polynomial of degree at most 
two, then u is actually equal to u £ for all e. The global well-posedness for (|1.17|) 
has been established in [7j. We first prove that u yields a good approximation for 
u e on bounded time intervals: 

Proposition 1.7. Let d > 1, a > cr < 2/(d - 2) if d > 3, and a S 5(R d ). 
Lei u G C(R;S) oe i/ie solution to (|1.17p . and /e£ 

(1.18) ¥» e (t, z) = ^ d/4 n (t, X ^ (f) ^) e *W)+«tM»-*(t»)A 
Lor all T > (independent of e > 0), we have 



If in addition a > 1/2, 



sup U £ (t) - </(i)IUw) =0(VZ). 

0<t<T 



sup ||V s (t)-V s (*)lh =0(Vi). 

Remark 1.8. The presence of it, which solves a nonlinear equation, clearly shows 
that the nonlinearity modifies the coherent state at leading order. Note however 
that the Wigner measure of i/j £ (see e.g. [191132] ) is not affected by the nonlinearity: 

w{t,x,0 = \\u(t)\\ 2 L2(Ild) 6(x-x(t))®5(Z-Z{t)) 

= \\af L2(ni) S(x-x(t))®5(Z-at))- 
The Wigner measure remains the same because the nonlinearity alters only the 
envelope of the coherent state, not its center in phase space. 
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Remark 1.9 (Supercritical case). Consider the case a < a c , and assume for instance 
V = 0. Resuming the scaling (|1.8|) . Equation (|1.15|) becomes 

id t u e + ^Au e = \e a ~ a "\u e \ 2,T u £ . 

At time t = 0, u € is independent of e: u^ t=Q — a. Setting ti 2 — e a "~ a and changing 
the time variable to s = t/H, the problem reads 

n 2 

(1.19) iHd s u H + — Au H = \u h \ 2<T u h ; u H (0, x) = a(x). 

Therefore, to understand the asymptotic behavior of it as e — > (or equivalently, as 
h — > 0) for t G [0, T], we need to understand the large time (s G [0, T/hj ) behavior 
in (|1.19p . This corresponds to a large time semi-classical limit in the (supercritical) 
WKB regime. Describing this behavior is extremely delicate, and still an open 
problem; see [5J. 

In order to prove the validity of the approximation on large time intervals, we 
introduce the following notion: 

Definition 1.10. Let u G C(R;£) be a solution to (|1.17|) . and k G N. We say 

that (Exp)k is satisfied if there exists C = C(k) such that 

Va,/?GN d , |a| + |/3|<fc, \\x a d%u(t)\\ L2{Rd) < e ct . 

Note that reasonably, to establish (Exp)k, the larger the k, the smoother the 
nonlincarity z i— > \z\ 2a z has to be. For simplicity, we shall now assume a G N. 

Proposition 1.11 (From [7]). Let d < 3, a £ N with a = 1 if d = 3, a G 5(R d ) 

and G N. T/jen (Exp)k is satisfied (at least) in the following cases: 

• a = d = 1 and A G R (cubic one- dimensional case). 

• cr 2/d, A>0 and Q(t) is diagonal with eigenvalues u>j(t) ^ 0. 

• cr 2/d, A>0 and Q(i) is compactly supported. 

It is very likely that this result remains valid under more general assumptions 
(see in particular [3 §6.2] for the case a = 2/d). Yet, we have not been able to prove 
it. Let us comment a bit on these three cases. The first case is the most general one 
concerning the potential V and the classical trajectory x(t): the only assumption 
carries over the nonlinearity (the important aspect is that it is L 2 -subcritical). The 
other two cases concern L 2 -critical or supercritical defocusing nonlincarities. In 
the second case, V is required to be concave (along the classical trajectory), and 
the last case corresponds for instance to a compactly supported HessV, when the 
classical trajectory is not trapped. In this last case, we have actually better than an 
exponential decay: Sobolev norms are bounded, and momenta grow algebraically 
in time. The following result could be improved in this case. 

Theorem 1.12. Let a G 6>(R d ). If (Exp)^ is satisfied, then there exist C, C2 > 
independent of e, and eo > such that for all e G]0, £0], 

\\i> s {t)-ip s {t)\\ n < Viexp(exp(C 2 *)), ^ t < Cloglog ~ 

e 

In particular, there exists c > independent of e such that 

sup \\^(t)-<p e (t)\\ n — >0. 

O^t^cloglogi 
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In the one-dimensional cubic case, this result can be improved on two aspects. 
First, we can prove a long time asymptotics in L 2 provided (Exp)z is satisfied. 
More important is the fact that we obtain an asymptotics up to an Ehrenfest time: 

Theorem 1.13. Assume d = er = 1, and let a S <S(R). If (Exp)?, is satisfied, then 
there exist C, C3 > independent of e, and £0 > such that for all e €]0, £0], 

||V £ (i)-V> £ (i)IU*(R) <x/iexp(C 3 i), 0<i<Clogi. 
In particular, there exists c > independent of e such that 

sup \\r(t)-^(t)\\ LHn) — ► O. 

If in addition (Exp)$ is satisfied, then for the same constants as above, 

U s {t) - V e {t)\\n < V^exp(C 3 t), O^t^ Clog-, 

e 

and 

sup \\r(t)-v s (t)\\H— >o. 

O^t^clogi £ ^° 

The technical reason which explains the differences between Theorem 11.121 and 
Theorem II .131 is that the one-dimensional cubic case is L 2 -subcritical. This aspect 
has several consequences regarding the Strichartz estimates we use in the course of 
the proof. 

These nonlinear results are to be compared with previous ones concerning the 
interaction between a linear dynamics (classical trajectories) and nonlinear effects. 
Consider the WKB regime 

(1.20) ied t ip E + —Aifj £ = V(x)i(; £ + \\ip e \ 2a i; £ ; ^ £ (0, x) = e^a{x)e ix ^°/ e , 
with V satisfying Assumption ll.il Like above, it is equivalent, up to a rescaling, to 

F 2 

ied t tjj E + —A^ £ = V{x)il> £ + \e a \^ £ \ 2(T ^ £ ; ^ £ (0,x) = a{x)e ix ^°/ e , 

with a = 2cr/3. The critical value in this regime is a c = 1 (see [6]). In (|1.20[) . 
this corresponds to initial data of order £ 1 /( 2cr ) in L°° , like in the present case of 
wave packets. However, the critical nonlinear effects are very different in the case 
of (fOQ|) . The following asymptotics holds in L 2 (R d ) (see [6]): 

ip £ (t,x) - a{t,x)e' 9(t ^e^ x ^ £ , 

e— >0 

as long as the phase </>, solution to the Hamilton-Jacobi equation 
d t <f>+^\V^\ 2 + V = ; <f>(0,x)=x-tD, 

remains smooth. More general initial phases are actually allowed: we consider an 
initial phase linear in x for the comparison with (|1.3j) . The amplitude a solves 
a linear transport equation: at leading order, nonlinear effects show up through 
the phase modulation g (which depends on A and a). This result calls for at least 
two comments. First, this nonlinear effect is rather weak: for instance, it does not 
affect the main quadratic observables at leading order, \tp £ \ 2 (position density) and 
elmip \7ip £ (current density). In the case of (|1.3j) . the profile equation is, in a sense, 
more nonlinear, even though in both cases, Wigner measures are not affected by 
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the critical nonlinearity. Second, the validity of WKB analysis is limited in general, 
even if V is a polynomial. If V = 0, <fi(t,x) — x ■ £o — i|£o| 2 /2 is smooth for all 
time, a(t, x) = ao(x — t£o) remains bounded, and the asymptotics can be justified 
up to Ehrcnfcst time, by simply resuming the proof given in [6]. If V(x) = E ■ x, 
Avron-Herbst formula shows that this case is essentially the same as V — 0. On 
the other hand, if V(x) = uj 2 \x\ 2 /2, classical trajectories in (|1.5j) are explicit: 

sm(u>t) 

x(t) = x cos(ui) + £ . 

CO 

They all meet at £o/ w at time <* = n/ (2u>): the phase <fi becomes singular as t — ► t*, 
and WKB analysis ceases to be valid, while the wave packets approach yields an 
exact result for all time in such a case. 

In [51 Q31 [T71 [2H HH1 > the authors have considered a similar problem, in a 
different regime though: 

(1.21) ied t ip s + y a *P £ = V(x)ip £ -\itj £ \ 2 < 7 ip e ; ip e (Q,x) = Q (^—y^j e i£o ' x/E , 

where Q is a ground state, solution to a nonlinear elliptic equation. They prove, 
with some precision depending on the papers: 

^(t,x)^ o Q ( j'^ ^ e^W'^W), 6 e (t) £ R. 

As pointed out in [24j , such results may be extended to an Ehrenfest time. An 
important difference with our paper must be emphasized, besides the scaling: the 
particular initial data makes it possible to rely on rigidity properties of the solitary 
waves, which do not hold for general profiles. In [S], some results concerning a 
defocusing equation with more general initial profiles are proved (or cited), in the 
same scaling as in (| 1 . 2 1 [) : however, it seems that unless V is a polynomial of degree 
at most two, only partial results are available then (that is, on relatively small time 
intervals). Finally, even when d 7 V = for all I7I ^ 3, the time intervals on which 
some asymptotic results are proved must be independent of e. 

1.3. Nonlinear superposition. We still suppose a — a c . For simplicity, in this 
paragraph, we assume that a is an integer: this is compatible with the fact that 
the nonlinearity is energy-subcritical only if d ^ 3. We consider initial data corre- 
sponding to the superposition of two wave packets: 

with <zi,a 2 S 5(R), (2:2,6) G R 2 , an d ^ (^2,6)- Fot 3 & I 1 , 2 }, 

(xj(t),£j(t)) are the classical trajectories solutions to (|1.5p with initial data (xj,£j). 
We denote by Sj the action associated with (xj(t),£j(t)) by ()1.7() and by Uj the 
solution of (|1.17| for the curve Xj(t) and with initial data aj. We consider <p e - 
associated by (|1.18[) with Uj,Xj,£j,Sj, and ip e £ C(R; E) solution to (|1.3[) with 
a = a c and the above initial data. 

The functional setting used to describe the function %p e must be changed in 
the case of two initial wave packets: recall that TC is defined through A £ and B £ , 
which are related to the Hamiltonian flow. The geometric meaning of A £ and B £ 
becomes irrelevant in the case of two wave packets. Instead, we use norms on S 
whose geometric meaning is weaker, since essentially, they reflect the fact that we 
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consider e-oscillatory functions, which remain somehow localized in space (before 
Ehrcnfcst time): 

ll/lk = ll/IUw) + WeVfWmm + IMU^)- 

For finite time, we have: 

Proposition 1.14. Let d < 3, a G N (a = 1 if d = 3), and a x ,a 2 G 5(R d ). For 
all T > (independent of e), we have, for all 7 < 1/2: 

sup W{t)-vi(t)-ipi(t)\\v e = o{e<), 

Besides, nonlinear superposition holds for large time (at least) in the one-dimensional 
case, if the points (£i,£i) an d (#2, £2) have different energies. 

Theorem 1.15. Assume that d = 1, a is an integer, and let ai,<22 G <S(R). 
Suppose that E\ ^ E2, where 

Suppose that (Exp)k is satisfied for some k 4 (for u\ and ui). 

1. There exist C, C3 > independent of e, and e$ > s«c/i £/ia£ /or all e G]0,£o], 

||V £ W - fi(tf - ^(t)||c. < £ 7 exp («p(C 3 t)) , < t < C log log J, 
wii/i 7 = 2 fc fc 7. 2 2 • iH particular, there exists c > independent of e such that 
sup ||^ c (t) - v?i(t) e - ¥^(t)||n B — *0. 

2. Suppose in addition that a = 1. There exist C, C4 > independent of e, and 
Eq > swc/i i/ia£ /or o/Z e G]0, £0], 

||^ e (t)-*»i(*) e -^(*)|| E . <e 7 e C4t , 0<t<Clogi «^ 7 =AzA. 
In particular, there exists c > independent of e such that 

sup ||V e (*) - Vi(*) e - — ; >o. 

It is interesting to see that even though the profiles are nonlinear, the super- 
position principle, which is a property of linear equations, still holds. There are 
many other such nonlinear superposition principles in the literature, and we cannot 
mention them all. 

This result is to be compared with those in [21] (see also references therein) , for 
several reasons. In [51] , the authors construct a solution for the three-dimensional 
Schrodinger-Poisson system which behaves, in 7f 1 (R 3 ) and asymptotically for large 
time, like the sum of two ground state solitary waves. The two solitary waves are 
centered, in the phase space, at the solution of a two-body problem: unlike what 
happens in our case, there exists an interaction between the trajectories, due to the 
fact that the Poisson potential is long range. In our case, the long range aspect 
of the nonlinearity (when d = a = 1; see [M]) does not have such a consequence: 
we will see that the key point in the proof of the above two results is the fact that 
in the wave packet scaling, the two functions ip\ and <p\ do not interact at leading 
order in the limit e — *■ 0: the nonlinear effects concentrate on the profiles, along 
the classical trajectories, and it turns out that these trajectories do not meet "too 
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much". In |8J, another nonlinear superposition principle was proved, in the scaling 
of (ll.2ip . However, in 8 , nonlinear effects were localized in space and time, so 
most of the time, the nonlinear superposition was actually a linear one. 

1.4. Outline of the paper. In Scction[2l we first analyze the linearizable case and 
prove Proposition 11.51 after a short analysis of the linear case. Then, in Section [3j 
we recall basic facts about Strichartz estimates in this semi-classical framework 
and prove the consistency of our approximation on bounded time intervals. Theo- 
rem ll.l2l is proved in Section[3J Finally, Section[5]is focused on the one-dimensional 
cubic case and Section [6] on the analysis of the nonlinear superposition. 

Notation. Throughout the paper, in the expression e ct , the constant C will denote 
a constant independent of t which may change from one line to the other. 

2. The linearizable case 

In this section, we assume a > a c and we prove Proposition II. 51 We first recall 
estimates in the linear case A = which are more precise than in £11.11 

2.1. The linear case. We suppose here A = 0. The first remark concerns the 
properties of the profile v. It is not difficult to prove the following proposition. 

Proposition 2.1. Let d ^ 1 and a £ S(R d ). For all k £ N, there exists C > 
such that the solution v to (jl.lip satisfies 

Va,/3eN rf , |a| + |/J|<ft, ||^^(i)|| L2(Rd) < e ct . 

A general proof of Proposition 12. II is given for instance in §6.1]. Let us now 
consider w € = ip £ — ipf in . We have u> e (0) = and 

ied t w e + -Aw e = V{x)w e - {V(x) - T 2 (x,x{t))) ^f in , 

where T 2 corresponds to a second order Taylor approximation: 

T 2 (x,a) := V(a) + (WV(a), x - a) + - (Hcss^(a)O - a), x - a) . 

We have seen in ^l.ll that the standard L? estimate for Schrodinger equations yields 

|K(*)|U2 (RIi) < V^\\y 3 v{t)\\ L 2 {Ild) < ^~ee ct . 
In order to analyze the convergence in E, we can write 

s 2 



isd t + —A - V(x) ) (eVw e ) = eWw £ - eVL e 



e 2 



: ,ed t + yA- V{x) ) (xw e ) 



yA,, 



xL £ = e 2 Vw £ - xL e , 



where 

(2.1) L £ {t, x) := (V(x) - T 2 (x, x(t) j) <^f in (i, x). 

Typically if d = 1, 

1 f 1 /// 

L £ {t, x) = -(x- x(t)) 3 ^(t, x) V (x(t) +6(x- x(t))) 9 2 dd 

_ (x~x(t)) 3 ^_ i(s(t ) + £ (t ).( x _ x(t)) y e _ {^ X — X(t) 



2eV4 
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where 

I(x,x(t)) = [ V" (x(t) + 9{x- x(t)))6 2 d6. 
Jo 

Energy estimates make it possible to show 

||eVu; £ (i)|| L 2 (R£i) + \\xw £ (t)\\ L 2 {Rd) < ^fee ct . 

However, the operators A £ and B £ defined in the introduction yield more precise 
results. For instance, || eV^iin || is of order 0(1) exactly, because of the phase 
factor in (|1.12[) . We note the formula 

(2.2) A £ (t) = ^iV-4^ = V-< 
\ e 



i{S{t)+i{t)-{x-x(t)))/e 



V e 



-i(S(t)+$(t)-( a: -x(t)))/ e . 



so for instance ||^4 e (£)<y2iin||i 2 is of order 0(1): morally, we have gained a factor %fe. 
Lemma 2.2. The operators A £ and B £ , defined by 

A*{t) = y/2V-M ; B £ {t) = X —j^, 



satisfy the commutation relations: 



ied t + -A-V,A s (t) 



ied t + -A~V : B £ (t) 



= ^i(VV(x)-VV(x(t))), 
= eA s (t). 



We can then write 

-2 



(ied t + Y A ~ = (W(k) - W (»(*))) w £ - A £ (t)L £ 



ied t + — A - V{x) ) B £ (t)w £ = eA £ {t)w £ - B £ {t)L £ 



In view of (|2.2p . we observe 

\\A £ (t)L £ \\ L2[Rd) <e^ 2 \\x 2 v(t)\\ LHRd) + e 3 / 2 \\x 3 Vv(t)\\ LHRd) + e 2 \\x 3 v(t)\\ LHRd) 

< e 3 l 2 e c \ 
thanks to Lemma [2TT1 Similarly, 

\\B £ {t)L £ \\ LHRd) <e 3 ' 2 e ct . 
Since we have the pointwise estimate 

| Ve (V7(:c) - W (x(t)))w £ \ < e \B £ (t)w £ \ , 
energy estimates yield 

\\™ e (t)\\H< [ {\\w £ (s)\\ H + V^e Cs )ds. 



We conclude by Gronwall Lemma: 

We will see in the following subsection that the arguments are somehow more 
complicated in the nonlinear setting. 
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2.2. Proof of Proposition 11.51 We now assume A ^ 0, and a > a c . For the 

simplicity of the presentation, we give the detailed proof in the case d — 1 only. 
We set again w e = ip e — ipf in and we write the equation satisfied by w s : 

ied t w e + ^d 2 x w £ = V(x)w s -(V(x)-T 2 {x,x{t)))ipf in + N s ; «,f t=0 = 0, 
where the nonlinear source term is given by 

First, since Xe a \ipf in + w £ \ 2a S R, the L 2 energy estimate for w £ yields 

l|w £ (£)|| jL 2 (R . ) < -||i E ||z,i([0,i];I,= (R)) + - W^Min + ^ I^Vfin || L i ([0,t];^(R)) ' 

where we have kept the notation (|2.ip . The contribution of N e cannot be studied 
directly, since we do not know yet how to estimate w e : since w e will turn out to be 
small, we use a bootstrap argument. 

Since we have llvfmWIIi^tR.) = e_1 ^ 4 || w (*)IU 00 (R)> Proposition ^. II and Sobolev 
embedding show that there exists Co > such that 

||vf in (t)|U«. (R ) ^ C ^ 1/4 e Cot , Vi > 0. 

The bootstrap argument goes as follows. We suppose that for t € [0, r] we have 

(2.3) IKWIU- ^ £- 1/4 e Cot , 

with the same constant Co. Since U>f t=0 = and i/) s €E C(R; S), there exists r e > 
(a priori depending on e) such that (|2.3p holds on [0,r e ]. So long as (|2.3[) holds, 

F IVlin+W I Vlin|| L 2 (R) < £ 7 ||a||i»(R)e . 

We infer 

lk e (t)IU'(R) < Vie ct +e Q -^ e 2CTCot . 
Applying the operators A e and £> e to the equation satisfied by w e , we find: 

f ie<9 t + y a 2 _ v(x)j A £ w £ = y/i (V'(x) - V (x(t))) w £ - A e L e + A £ N e , 

e 2 \ 
\ed t + —d 2 x - V{x) J B e w e = eA e w e - B e L £ + B £ N e . 



We observe that in view of (12. 2|) . A £ acts on gauge invariant non linearities like a 
derivative. Therefore, so long as (|2.3|) holds, 

\\A s (t)N £ (t)\\ L2{R) < e<* (||^n(*)lli»(R) + IK(i)llxr-<H)) ll^(*)vU*)IU»<H) 

+ ^ (lltfi„(t)ll£~(R) + IK(*)ll!~(R)) ||A £ (t)^(t)|| i2(R) 
< £ a- CT /2 e 2 CT Cot ( e Ct + ||^ Wu;£W || l2(r) ) _ 

Similarly, we obtain 

\\B%t)N%t)\\ LHn) < e-^e 2 ^* (e CT + (t)|| za(R) ) . 
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We infer, thanks to the linear estimates, 

\\A e (t)w e (t)\\ L 2 (n) < \\B e w £ \\ L i m] , L 2 {Il)) +^e ct 

(e Cs + \\A%s)w% S )\\ LHR) )d S , 



t 



\B £ (t)w s (t)\\ L 2 {R} < \\A £ w £ \\ LHm , L 2 {R)) + ^fee ct 

3 2aC s ( e Cs + llB e is)w s {s)h2{R ^ d ^ 

Gronwall Lemma yields, so long as (|2.3[) holds: 

\\w £ {t)\\n < J £ 7 e Cs exp ^Ce a ~ a " J e 2aC ° s ' ds' \ ds 

< exp (Ce a ~^e 2aCot ) f e^e Cs ds < exp (Ce a ~ a ^e 2aCot ) e 7 e ct , 
Jo 

where 7 = min(l/2, a — a c ). First, we notice that 

£ a- ace 2*C t <± for ^ ^ lQg I. 

2aC e 

Then, setting k — f^ffi , Gagliardo-Nirenberg inequality yields, so long as 
holds, with also t ^ Klog ~, and thanks to the factorization (|2.2j) , 

IK(*)IU-(R) Z ^II^WIIl^II^W^WII^r) <^ 1/4 e ct . 

This is enough to show that the bootstrap argument (|2.3[) works provided the time 
variable is restricted to 

Ce^e ct < 1, 

that is, ^ t ^ Clogi for some C > independent of e. Proposition 11.51 follows 
in the case d = 1. 

To prove Proposition 11.51 when d 2, one can use Strichartz estimates. This 
approach is more technical. Since the case a > a c does not seem the most interest- 
ing one, and since we will use Strichartz estimates in the fully nonlinear case, we 
choose not to present the proof of Proposition 11.51 when d 2. 

3. Fully nonlinear case: bounded time intervals 



In this section, we prove Proposition 11.71 This gives us the opportunity to 
introduce some technical tools which will be used to study large time asymptotics. 

3.1. Strichartz estimates. 

Definition 3.1. A pair (q,r) is admissible if 2 ^ r ^ ^4; (resp. 2 ^ r ^ 00 if 
d = \, 2 ^ r < 00 if d = 2) and 

q { ' \2 r 

Following [20] 140] [27] . Strichartz estimates are available for the Schrodinger 
equation without external potential. Thanks to the construction of the parametrix 
performed in |15[ I16j , similar results are available in the presence of an external 
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satisfying Assumption 11.11 (V could even depend on time). Denote by U £ (t) the 
semi-group associated to — ^-A + V: 4> £ (t, x) — U £ (t)(j)o(x) if 

Kftf + yAf = 7f ; ^ e (O,aj)=0o(a?). 

From [15] , it satisfies the following properties: 

• The map 1 1— » U £ (t) is strongly continuous. 

• U £ (t)U £ (s) = U £ (t + s). 

• u E (ty = u e {t)- 1 = u e {-t). 

• U E (t) is unitary on L 2 : \\U £ (t)(t>\\ L 2 (Rd) = \\<j)\\ L 2 iRd) . 

• Dispersive properties: there exist S, C > independent of e G]0, 1] such 
that for all t G R with \t\ < S, 

C 

II^ £ WIIli(r^l~(r^) «S ( £ \t\)d/2 - 
We infer the following result, from [27] : 

Lemma 3.2 (Scaled Strichartz inequalities). Let (q,r), (qi,ri) and (92,^2) be ad- 
missible pairs. Let I be some finite time interval. 

1. There exists C — C(r, \I\) independent of e, such that for all <j> £ L 2 (R. d ), 
(3.1) e 1 ^ ||I^(-)0ll £ « ( j.£r (Rd) ) < CU\\ L2{Rd) . 

2. Lf I contains the origin, S I, denote 



D s I (F)(t,x) = U £ (t- s)F(s,x)ds. 

Jln{s^t} 

There exists C = C(ri,r2, |/|) independent of e such that for all F £ L q2 (I; I/ 2 ), 

(3.2) e V*+i/» Pf(F)|| L51(/;i . 1(Rd)) <C||F||^ (z;i4(Ri)) . 

3.2. Proof of Proposition IT771 Denote the error term by w £ = -0 e — <fi £ , where 
ip £ is now given by (|1.18p . and u £ C(R; S) satisfies (|1.17j) . This remainder solves 

(3.3) ied t w £ + —Aw £ = Vw £ - C £ + \e a " (\ip £ \ 2a ip £ - |<^| 2< V) ; wf t=0 = 0, 
where 

(3.4) C £ (t,x) = {V(x)-T 2 (x,x{t)))<p £ (t,x) 

is the nonlinear analogue of L £ given by (|2.ip . Duhamel's formula for w £ reads 

ft + T 

w £ (t + t) = U £ {T)w £ {t) + ie^ 1 J U £ (t + r- s)C £ (s)ds 

ft + T 

- iXe "- 1 / U £ (t + t-s) {W\ 2a ^ £ - \p £ \ 2a <p e ) {s)ds. 



Introduce the following Lebesgue exponents: 

_ 2ct(2ct + 2) 4ct + 4 



2 - (d - 2)ct ' * do- 
Then (g, r) is admissible, and 

1 2cr 1 12crl 



r = 2o + 2. 
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Let t 5* 0, r > and I = [t, t + r]. Lemma f3 . 2 1 yields 

+ e^- 1 - 2 /" \\\r\ 2a r - l^| 2 > e |L 9 ' (/; ^) 
In view of the pointwise estimate 



we infer 



(3.5) 



£ a c -l-2/q 



Thanks to [7], we know that the rescaled functions for ip £ and ip £ , are such that 
u e , u € C(R;£), with estimates which are uniform in e g]0, 1]. Typically, for all 
T > 0, there exists C(T) independent of e such that 

ll^« e |U->([0,n;£») + II^IU~([0,T];£*) < C(T), P £ {Id, V, x}. 

In terms of i\) e and ip s , this yields 

(3.6) ||7> £ V1l~([o,t];L 2 ) + ||P> £ |U-([o,r];^) < C(T), V e € {Id, A £ , 5 £ }. 
The formula (|2.2[) and Gagliardo-Nirenberg inequality yield, if ^ S(p) < 1, 

(3.7) \\f\\ LP(Rd) ^^L ll/ll^g m e (*)/ll^), V/ei? 1 ^), VteR. 
We infer that there exists C(T) independent of e such that 

(3.8) ||^ £ WIU-(R d ) + ||^ £ WI|l,(r,) < C(T)e- s ^/ 2 , Vi € [0,T]. 
Recalling that / = [t, t + r], we deduce from (|3.5p : 

+ e a e -l- a/9ra <r/ fle - rfW | K || t<(j!ir)< 

Since (q, r) is admissible, we compute 

3.9 a c - 1 rf r = -= = 0, 

hence 

(3 ' 10) +^11^11^. 

Choosing r sufficiently small, and repeating this manipulation a finite number of 
times to cover the time interval [0, T], we obtain 

(3.11) IKH^acT];^) < e" 1/9 ||^ £ |Ui([0,T];^) + <r 1_1/g II C e || L i ([0 , T ] ; L>) ■ 

Using Strichartz estimates again and (|3.8j) . we have, with J = [0, i] and ^ t ^ T, 

< £ - 1 ||/: £ ||l 1 (,/;l 2 ) + ^- 1 - 1/9 ||I^ £ I 2ct V' £ -|^ £ I 2 V £ ||^ (J; ^ ) 
^ I|w £ ||li(,/ ; l2) +e _1 ||/: £ || L i(./ ; L2) 

I _a c — 1— 1/g 1— 1/g— er<5(r) || /•£ | 

^ I|w £ ||li(,/ ; l2) +e~ 1 ||£ £ ||li(./;L2), 
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where the last estimate stems from (|3.9[) . We have the pointwise control 
\C £ \ <\x- x(t)\ s \<p £ (t, *)\ = e 3/2 e~ d/i (\y\ 3 \u(t, y)\) 

We infer 



■w 



7? 



-In re 



-C e ||Ll([0,T];L2(Rd)) ^V^IIM u(i,y)|| L l([ ,T];L2(R<*)), 



and the first part of Proposition 1 1 . 71 follows from Gronwall lemma. 

To establish a control of the 7i-norm, we notice that in view of (12. 2 [I . we have, 
for 7> e G 

7> £ (|<£ £ | 2<T <£ £ ) pa |<f | 2CT 7? £ <f , 

where the symbol "ps" is here to recall the abuse of notation when P £ = yl £ (there 
should be two terms on the right hand side, with coefficients). Lemma [2.21 shows 
that we have 

(isdt + yA-yj A £ w £ = y/i (VF(x) - W (x(t))) w £ - A £ C £ 

+ \e a °A e {m 2 °r-\v e ?> e )- 

The first term of the right hand side is controlled pointwise by Ce\B £ w e \. The 
L 2 -norm of the second term is estimated by 

\\A £ (t)c £ (t)\\ L2(Rd) < e 3 / 2 (\\\y\Mt,y)\\L^) + IIM 3 v«(i, y )|| i2(Rrf) ) . 

Finally, we have 

A £ (\^ e \ 2a ip e - \<p s \ 2a v s ) « W\ 2a A £ ij £ - \ v e \ 2a A e <p e 

m \w £ + ip £ \ 2a {A £ w £ + A £ tp £ ) - \ip £ \ 2a A £ Lp £ 
(3.12) « K + <p £ \ 2a A £ w £ + (\w £ + (p £ \ 2a - \(p £ \ 2a ) A e ip £ . 

The first term of (|3.12[) is handled like in the first step. For the second term, we 
have, since a > 1/2, 

\\w £ + V > £ \ 2 < 7 -\^ £ \ 2a \ < (\vf\ 2 °- 1 + \<ff\* r - 1 )\vf\. 

Following the same lines as for the L 2 estimate, we find 

\\A £ w £ \\ Lq(I , Lr) < e-^WA-^it)^ + e- 1 ^\\B £ w £ \\ L1(I , L , ) 



+e~ 1 ~ 1 /i\\A £ £ £ \\ Ll(IlL2) + t 2 °' s \\A £ w £ \\ Lq{hLr) + r 2 °' s 



L*(I;Lr) T ' 11^ \\Ll{I;Lr) i 



by using the estimate 



\\A*{t)<p*mL*M<e- S{r)fi \^ 



and the remark 

||A e (i)V(t)IU»{R<) = l|V 2 M(<)|| L 2 (R£i) . 

Since cr > 1/2, the nonlinearity z i— > |z| 2cr z is twice differentiable, and one can prove 
u e C(R; H 2 (R d )) (7). Using (|3.11|) and the same argument as in the first step, 
we infer 

\\A £ w £ \\ Lq{I . iLn < e- 1 ^\\A e (t)vf(t)\\L' + £- 1/<z ||i?W|| L1(/;L2) 

+ E-^v^w^,^ + s- 1 - 1 '* ||/: £ || L1(/;i2) , 
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hence, using Strichartz estimates again, 

WA^Uoo^ < \\A £ (t)w s (t)\\ L 2 + ||BW|| L1(/;L2) + e^WA^Wli^ 
+ e^\\A*w^ L , iI]Lr) +e^\\w^ L , iI]Lr) 
< \\A e (t)w e (t)\\ L2 + ||B e to B || il(/ . £a) + e- 1 \\A e £?\\ LHl . iL , ) 



1/9 t 



It) 



lL9(/;f) 



Since we have similar estimates for B e w e , we end up with 

P E w E |U°°(J;£=) + ||B e «; e || £ « C j; La) < ||A £ u; e || L i (J . L 2 ) + ||B% e ||z,i(J;i2) 

< \\A £ w £ \\ L i (J , L 2 ) + + 
Proposition 1 1 . 71 then follows from Gronwall lemma. 

4. Fully nonlinear case: proof of Theorem 11.121 

To prove Theorem 11.121 the strategy consists in examining more carefully the 
dependence of the L Z7-norms with respect to time in the previous proof. Also, 
since (Exp) 4 concerns only u, not u E , we need a bootstrap argument in order to 
use the same control for the error term w E as for the approximate solution ip e . This 
control carries on the L r (R d )-norms, for fixed t. By (Exp)\, the relation 

\\A £ {t)<p £ {t)\\ L 2 (Rd) = \\Vu(t)\\ L 2 {Rd) , 

and the modified Gagliardo-Nirenberg inequality (|3.7D , we have the following esti- 
mate, for all time: 

(4-1) M(t)\\ L ^)<e- S{r)/2 e Kt . 

We will use the following bootstrap argument: 

(4.2) IK(i)IU-(R<) < e-'W/'c"*, t€[0,n 

with the same constant k as in (|4.1[) to fix the ideas. By Proposition 1 1 . 71 for any 
T > independent of e, (|4.2p is satisfied provided < e ^ e(T). By this argument 
only, it may very well happen that e(T) — > as T — > +oo. The goal of the bootstrap 
argument is to show that we can take T e = C log log \ for some C > independent 
of e, provided that e is sufficiently small. 

The key step to analyze is the absorption argument, which made it possible to 
infer p. lip from ()3.10|) . We resume the computations of £13.21 from the estimate 
(|3.5[) . Rewrite this estimate with J = [t, t + r], t, r ^ 0: 

lk £ IU«(w) ;$ ^ 1/9 lk £ (*)IU= + e-'-^wc^^j,^ 

+ e a "- 1 - 2/ ' 1 (\\w%\ I;Ln + ll^lllV^)) IKIIt^) • 

For simplicity, assume r ^ 1: (|4. 1[> and (|4.2|) yield, in view of (|3.9|) . 

||w £ || L ,(/^) < M(e- 1 /«K(i)|U 3 + £ - 1 - 1 /«|| J C £ |U 1(J;i2) 
+ r i/« e 2- t || w s|| iq(/;ii . 
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for some constant M independent of e, t ^ and ^ t ^ 1. In order for the last 
term to be absorbed by the left hand side, we have to assume 

Mr^e"^, that is, t ^ Ce- Ct 

for some C independent of e, t ^ and ^ r ^ 1. Proceeding with the same 
argument as in §3.2[ we come up with: 

< / P Cs \\v, e \\r^nn .l.r^ds + F- 1 [ P Cs \ 



\\w £ \\ L ^([o,t];L--) < / e Os ||w £ || L oo ([0 . s];L 2 ) ds + £ / e° s ||£ e (s)|| L 2ds 
Jo Jo 

< f e Cs \\w £ \\ L ^ {[0 . s] . L 2 ) ds + y^e c \ 
Jo 

where we have used (Exp)^. Gronwall lemma yields: 

\\w £ \\l^([o, t];L2) ^ \/eexp (Cexp(Ct)) < y/eexp (exp(2Ci)) . 
Mimicking the computations of £13.21 we have, thanks to (Exp) 4 and so long as (|4.2 



holds, 

\\A e w e \\ L ^ { [o,t};L^ + \\B £ w £ \\ L ^ {[Q ^ t] . L 2 ) ^ \/eexp (exp(C<)) . 

To conclude, we check that (|4.2p holds for t ^ c log log |, provided c is sufficiently 
small. Gagliardo-Nirenberg inequality (|3.7p yields 

ll^ e (0llz.^C^) < ^-^^^ll^ll^ro^l^^ll^^lll^ao,*];^) 

M£- 5 M /2 Veexp (exp(Ct)) . 

Therefore, taking e sufficiently small, (|4.2[) holds as long as 

M^exp(exp(Ct)) «S e Kt . 

We check that for large t and sufficiently small e, this remains true for i ^ clog log i, 
with c possibly small, but independent of e g]0,£o]- This completes the proof of 
Theorem [HI 

5. Ehrenfest time in the one-dimensional cubic case 

As pointed out in the introduction, since we consider nonlincaritics of the form 
z 1 — > \z\ 2a z with cr e N, the one-dimensional cubic case is special. Not because 
it is integrable (see Remark 11.61 (|1.17p is not completely integrable, unless no 
approximation is needed to describe the wave packets, ip £ = tp £ ), but because it 
is the only case where the nonlinearity is L 2 -subcritical, a < 2/d. This case is 
in contrast with the general case of energy-subcritical nonlinearities: without any 
other assumption on Q(t) than Q £ C°°(R; R) n L°°(R), it seems that the only a 
priori control that we have for u, solution to (jTTTTJ) , is 

(5.1) ||u(f)|| L 2 (R d) = ||a|| L 2 (R d), Vi £ R. 

A remarkable case where other a priori estimates are available is when Q is constant, 
but in this case, ijj £ = Lp £ . Otherwise, the most general reasonable assumption seems 
to be (Exp)k, which has been considered in the previous section. Note also that if 
d = 1, the notations of £13.21 become: 



3 , * ~ , - 

So to improve the result of Theorem 1 1.1 21 we assume a = d = 1 and start with the 
crucial remark: 
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Lemma 5.1. Suppose a = d = 1, and for a € L 2 (R), consider u G C(R; L 2 (R)) 
the solution to (|1.17|) . Then there exists C such that 

ll M llL s ([t,t+l]:L 4 (R)) < C||a|| L 2( R ), \/t E R. 

Proof. First, recall that since a = d = 1 and a € L 2 (R), (|1.17|) has a unique 
solution 

«eC(R;I 2 (R))nlf oc (R;i 4 (R)). 
In addition, (|5.ip holds. Denoting 

W(t,x) = ^V" (x{t))x\ 

it has been established in [7] that since V" € L°° (R; R) , uniform local Strichartz 
estimates are available for the linear propagator. Following [HI [16], let U(t,s) be 
such that as u(t,x) = U(t, s)uo(x) is the solution to 

idtu + — Au = W(t, x)u ; u(s, x) = Uo{x). 

Then Lemma l3~2l remains true (with e = 1) when U e (t— s) is replaced with U(t, s), 

t, s e r. 

Let t, t ^ 0, with r ^ 1, and denote / = [t, t + r]. Strichartz inequalities yield: 
||""]Us(J;i*) «S C(r)||u(i)|| L 2 + C(r) || \u\ 2 u\\ l8/7{i , l4/3) . 

In view of (|5.ip . and using Holder inequality after the decomposition 

3 3 1 7 3 1 

4 ^ 4 + ' 8 ~ 8 + 2' 

we infer 

< C(T)\\a\\ L 2 + C{T)y/¥\\u\\ 3 L B {I . L *y 

Since r 1, we may assume that C(r) does not depend on t: 

I|w||l8( /;L 4) C||a|| L 2 + CVr||tt|||,8(7 ;i 4). 

We use the following standard bootstrap argument, borrowed from J3J: 

Lemma 5.2 (Bootstrap argument). Let f — f(t) be a nonnegative continuous 
function on [0, T] such that, for every t € [0,T] , 

f(t)^M + Sf(t) e , 

where M, 5 > and 9 > 1 are constants such that 

M< { 1 ~e) (0«s)i/(«-i) 5 /(0) ^ (^)V(fl-i) • 

Then, for every t£ [0,T], we /icwe 

/(*) < 

Lemma IQl follows with [i, t + 1] replaced with [£, t + r] for < r ^ to <C 1. We 
then cover any interval of the form [t, t + 1] by a finite number of intervals of length 
at most tq, and Lemma 15. II is proved. □ 
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Proof of Theorem \1.13[ Like in the previous section, we resume the proof of Propo- 
sition [TT7J and pay a more precise attention to the dependence of various constants 
upon time. We modify the bootstrap argument of SJU in view of Lemma |5. 11 
is replaced by 

(5-2) IKIU«([M+i]iL«(R)) <£" V8 ||a|U=(R), Vie[0,T]. 



By Proposition 11.71 for any T > independent of e, (|5.2|) remains true provided 
< e ^ £(T). Keeping the notations of H3.21 we have: 

0=1 ; q = S ■ r = 4. 
With £ ^ 0, r e]0, 1] and I=[t,t + r], $M becomes 
\\w s \\l H i-M) < e- 1/8 |K(i)|U 2 + e-^WJ^Un!.^) 

+ £ V4 (11^1^8/3(^4) + ||<P S Hi8/3 (J;Z ,4)) IKIIz,8(I;/>) 

^e-^H^WIU.+e- 1 - 1 / 8 !!/:,!!^^) 

+ £l/4rl/4 (11^11^(7^4) + ||<y5 £ |liB (/ . L 4)) ||w e || Z .(Jj£* ) 

(5-3) < e-^II^WH^ + e-^^H^Hzi^) + r 1 ^ ||^|| L8(W) , 

where we have used Lemma 15.11 and (|5.2I) . Choosing r sufficiently small and inde- 
pendent of t, we come up with 

ll^ e |U°°([0,t];L2) < ||w e || L l ([0 ^ ];L 2) +£~ 1 ||£ £ ||Li([0,t];L2) 

< T-.irrn f.i-r.21 + \fe ( e Cs ds, 



\W \\m[0,t];L*) + V£ 

Jo 

by (Exp) 3. Gronwall lemma yields 

||t« S |U«([0,t];L») < V^e Ct . 

Back to (|5.3[) . we infer, with t<1, 

IKIU*(/ ;i *)<£ 1/4 e Ct . 

Therefore, there exists c > such that (|5.2p holds for T — c log ^ provided e is 
sufficiently small, hence the first part of Theorem 11.131 

It is then quite straightforward to infer the estimates in Ti., by rewriting the end 
of the proof of Proposition 11.71 with (|5.2[) in mind. □ 

6. Nonlinear superposition 

6.1. General considerations. The proof of Proposition ll.141 and Theorem 11.151 
follows the same lines as the proof of Proposition 1 1 . 71 and Theorem 1 1.1 31 The main 
difference comes from the way one deals with the nonlinearity, since new terms 
appear. These terms come from the nonlinear interaction between the two profiles 
ifi and if2- Denote w e — ?p e — <p\ — tp\. It solves 

£ 2 
2 

where we have now 



ied t w E + —Aw E = Vw e - C e + W £ ; iof t=0 = 0, 



£ E (t,x) = (V(x) -T 2 (x,x(t))) {<pl{t,x) + <pl(t,x)) , 



NONLINEAR COHERENT STATES AND EHRENFEST TIME 



21 



and 

Decompose J\f £ as the sum of a semilinear term and an interaction source term: 
J\f e = J\f§ + Aff, where 

Af§ = e<*° (\w £ + v £ 2 \ 2 °{w £ + + rt) - \cpl + V \\ 2 °(ip\ + tpl)) , 
Aff = e a < (\<pt + ri\ 2 °(ri + tf) - |^| 2 >! - |^| 2 VI) ■ 

We see that the term Afg is the exact analogue of the nonlinear term in (|3 . 3[) . 
where we have simply replaced <p 6 with Lp\ + ip\. We can thus repeat the proofs of 
Proposition 1 1 . 71 and Theorem ll.l3( respectively, up to the control of the new source 
term Aff (the linear source term C £ is treated as before) . More precisely, we have 
to estimate 

-\\Nf\\ L i(io, t ]-,L2(-R*))- 

The hrst remark consists in noticing that if a is an integer, Aff can be estimated 
(pointwise) by a sum of terms of the form 

To be more precise, we have the control, for fixed time, 

j\\Aff{t)\\ LHRi) < e da ' 2 J2 IMI' 1 x I^I^L'cr-) • 

We will see below why the right hand side must be expected to be small, when 
integrated with respect to time. We need to estimate 

xi(t) - x 2 {ty 



r (2(7/2 



L 2 (R d ) 



t, X 



u 2 2 (i, x) 



L 2 (R d ) 



with l\, £2 ^ 1, £\ + £2 = 2er + 1. We have the following lemma: 

Lemma 6.1. Suppose d ^ 3, and a is an integer. Let T £ R. < 7 < 1/2, and 

(6.1) i s (T) = {te [o,T], Mt)-x 2 (t)K^}. 

Then, for all k > d/2, 

\ £ \Wf(t)U e dt < (M fe+2 (T)) 2CT+1 (te^ 1 / 2 ^ + \I £ (T)\) e CT , 

where M k (T) = sup{|| (x) a d£ %-|| L ~([o,T];L 2 (R d ); 3 e {1,2}, \a\ + \[3\ ^ k} . 
Proof. We observe that for rj G R d , 



sup 

x£R d 



With rj e (t) = we infer (forgetting the sum over £i,£ 2 ), 



[0,T]\/»(T) 



L 2 (R 



dt < 



< 



< 



[0,T]\/ e (T) 



(x - r,*(t))- k (x)- k (x Tf (t))- fe t£ (t, x - rf(t)) u*"(t, x) 

f dt 

L°°(10,T]-,L*) J\0,T\\I*(T) \V £ (t)\ k ' 



L 2 



{x} U-ff 




/ \k l 2 




(x) u 2 2 




L=°([0,T]:L 4 ) 
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We have, for j £ {1,2}, 



/ \fe (3 
(x) u/ 



L°°([0,T];L i ) 



< 





k 




(x) 




Uj 




k 




(x) 




u 3 



>([0,T];L*) " 3 "L<~([O.T]xR*) 



L oo ([0,T];J? 1 ) 



l 3\\L°°([0,T];H k ) 



<M k+1 (T)^, 



where we have used i? 1 (R d ) C L 4 (R d ) since d ^ 3. On the other hand, 



< 



'[0,T]\I*(T) l?f 

On I e (T), we simply estimate 
1 



[0,T]\/-(T) -X 2 (t)| fe 



(ft < £ fc (V2-7) T . 



IK)' 11 ~ ll U l|ri°o ([ o, T ] xR )||M2||ioc( [0 ,T]xR)ll W 2||z,l(/ e m;i 2 (R)) 

<M fe (T) 2CT |/ e (T)||| U2 || L » ([0 , T];i2(R)) 

<M k {T) 2a+1 \r(T)\. 

The L 2 estimate follows, without exponentially growing factor. This factor appears 
when dealing with the E £ -norm. Typically, 

||£V^(t)|| i2(Rd) < v^llVuiCiJHi^a* 

< 



I (*) 1 1| «j (*) 1 1 i» (R") , 
H^iWIlL^R^) < <Je\\xuj{t)\\v>{n*) + \Xj(t)\\\Uj(t)\\ L z(B.<l). 
The result then follows from the above computations, and Lemma H~2l 



□ 



At this stage, the main difficulty is to estimate the length of I s (t). We do this 
in two cases: bounded t, and large time when d = 1. 



6.2. Nonlinear superposition in finite time. In the proof of Proposition 11.71 
we have only used the fact that u £ £ C(R; £), with estimates which are independent 
of e. Recall that in the case of a single wave packet, ip £ and u e are related through 
(ll.8[) : in the case of two wave packets, there is no such natural rescaling. So in 
the case of two initial wave packets, we are not able to prove uniform estimates 
for ip £ , like in (|3.6p . Even to prove Proposition 11.141 which is the analogue of 
Proposition ll.7l we need to use a bootstrap argument. We know that for j £ {1, 2}, 

H(*)IU-(r') < C(T) £ -*M/ 2 , vte[o,r]. 

The bootstrap argument is of the form: 

\\w s (t)\\ Lr(Ri) ^C(T)s- s ^ 2 , Vie [0,71, 

with the same constant C(T) if we wish. Repeating the computations of £)3.21 we 
first have, for t £ [Q,T] and so long as the above condition holds, 

l|w E |U°°([0,t];£ 2 ) < -H^IUhIO.T];!- 3 ) + - || Mj \\ L i ([0,T];L 2 ) ■ 

As we have seen in ij2.1[ (eWw £ ,xw e ) solves a system which is formally analogous 
to the system satisfied by (A £ w e , B e w £ ). Therefore, under the bootstrap condition, 
we come up with 



\\W ||L°°([0,t];£ £ 

We easily estimate 



<-H^I 



LHIO.T];^,) 



-ll-Yf ||L!([0,T];S e 
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so in view of Lemma 16. 1[ the point is to estimate the length of I e (T) . 
Lemma 6.2. For T > (independent ofe), we have 

\F(T)\=0(e"<), 

where I e (T) is defined by (16. ip . 

Proof. The key remark is that since ^ (2:2, £2)7 the trajectories x\(t) and 

X2(t) may cross only in isolated points: by uniqueness, if x±(t) — X2(t), then 
x\(t) 7^ X2(t). Therefore, there is only a finite numbers of such points in the 
interval [0,T]: 

(xi(-) -x 2 {-))- 1 (0)n [0,T] = {^Ji^j, where J = J(T). 

If we had J = 00, then by compactness of [0, T], a subsequence of (tj)j would 
converge to some r € [0, T], with x\(t) — X2(t). By uniqueness for the Hamiltonian 
flow, ±i(t) ^ X2{t): t cannot be the limit of times where x±(tj) = X2(tj). 

By uniqueness for the Hamiltonian flow, continuity and compactness, there exists 
S > such that 

J 

inf{|£i(i) - x 2 (t)\ ; t e 1(6, T)} = m > 0, where 1(6, T) = \J [tj - 6, tj + 6] , 

3=1 

and there exists e(6,T) > such that for e e]0, e(6, T)], I s (T) C 1(6, T). 
Let t e I s (T) n [tj - 6, tj + 6}. Taylor's formula yields 

xi(t) - x 2 (t) = xi(tj) - x 2 (tj) + (t- tj) (±i(t) - x 2 (t)) , t e [tj - 6,tj + 6}. 
We infer 

e 1 ^ \x x (t) - x 2 (t)\ > \t-tj\m, 
and Lemma [6721 follows. □ 

Back to the bootstrap argument, we infer 

ll^llL~([0, t ]; SE )<^+e fc(1/2 - 7) +^. 

Fix 7 e]0, l/2[. By taking k sufficiently large in Lemma [67T1 this yields 

l|w £ || L ~([o,t];S £ ) <e 7 - 
Gagliardo-Nirenberg inequality yields 

IKWIU" < E - 5W |k £ (t)||^ <5(r) || e V^(t)||i ( 3 r) ^e-^H^IU-dD.tjiS.) <e* r - tf(r) . 
To close the argument, we note 

£ 7-<5(r) e -*(r)/2 provided e < 1 and 7 > 

The last condition is equivalent to 7 > 4 ^ 4 , which is compatible with 7 < 1/2 
since the nonlinearity is energy-subcritical. 
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6.3. Nonlinear superposition for large time. Things become more compli- 
cated when T is large. We first need to control Mk' this is achieved assuming 
(Exp)k, and we have 

M k (t)<e ct . 

The main point is to estimate \I e \. This is achieved thanks to the following 
proposition, whose proof relies heavily on the fact that the space variable is one- 
dimensional. 

Proposition 6.3. Under the assumptions of Theorem \1.15\ there exist C,C'o > 
independent of e such that 

\I £ (t)\<e''e c ° t \E 1 -E 2 r 2 , 0<^Clog~ 

£ 

Proof of Theorem \1.15[ Before proving Proposition ^. 3[ we show why this is enough 
to infer Theorem 1 1.1 51 By Lemma 16. 11 we have, if (Exp)k is satisfied, 

Jl|JV/IUi([o, t ] ; s e ) <e Ct ( te ( fc - 2 )(V2-7) + £ Te Cot ) < ( e (*-2)U/2-7) + e ct_ 

Optimizing in 7, we require (k — 2) (1/2 — 7) = 7, that is 

k-2 
7 ~ 2k -2' 

We can thus resume the bootstrap arguments as in fjH and $5j respectively. The 
key is to notice that this works like in the previous paragraph, since 

k-2 a 
7= > (fc>4). 

This yields Theorem 1 1.1 51 □ 

Proof of Proposition \6.3[ We consider J E (t) an interval of maximal length included 
in I s (t) and N e (t) the number of such intervals. The result comes from the estimate 

\I £ (t)\ < N e (t) x max|J £ (i)|, 

with 

(6.2) \J s {t)\<e<e Ct \E l -E 2 \-\ 

(6.3) N £ {t) < te 2Ct \E x - E 2 \- x < e 3C, *|£i - E 2 \-\ 

We first prove prove (|6.2p . Let ti,t 2 € J £ (t). There exists t* € [iijia] such that 

\(xi(h) - x 2 {h)) - ( Xl (t 2 ) - X2 (t 2 ))\ = \t 2 - h\ lar ) - &(oi , 

whence 

\h-t 2 \^Mt*)~t 2 (t*)\- 1 x2e\ 

On the other hand, 

Mn - un\ > \Mn\ - \un\\ > ^^^pp - 

Using 

Mt*)\ 2 - |6(<*)| 2 = 2(E 1 -E 2 - V{ Xl {t*)) + V(x 2 (t*))) , 
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we get 

||6(**)| 2 -|6(**)| 2 |> \E t -E 2 \-eie c \ 

whence 

|*i -* 2 | <e^e ct \E 1 -E 2 \-\ 

provided e 1 e ct <C 1. 

Let us now prove ()6.3|) . We use that as t is large, N e {t) is comparable to the 
number of distinct intervals of maximal size where | x\ (t ) — x 2 (t ) | ^ e 7 . We consider 
J' £ = [t[ , t' 2 ] such an interval . We have 

Ixxit'J - x 2 (t' 1 )\ = \ Xl {t' 2 ) - x 2 (t' 2 )\ = e\ and We [4,4], \ Xl (t) - x 2 {t)\ ^ e\ 

Therefore, for t £ [t[ , t' 2 ] , the quantity x\ (t) — X2 (t ) has a constant sign: we suppose 
that xi(t) — x 2 (t) is positive. We then have 

6(4) - 6(4) > and a (4) - 6(4) < 0. 
Using the exponential control of V (x 3 -(£)) for j 6 {1, 2}, we obtain 

(6(4) - 6(4)) - (6(4) - 6(4)) < e ct |4 - 41- 

We write 

|I6(4)I 2 -I6(4)I 2 I 



6(4) -6(4) = 16(4) -6(4)1^ 



-6(4) + 6(4) = |£(* 2 )- 6(4)1 > 



16(4)1 + 16(4)1 

>e- Ct |l6(4)| 2 -|6(4)i S 

'I6(4)l 2 -I6(4)l 2 | 



16(4)1 + 16(4)1 



I6(4)I 2 -I6(4)P 



Besides, in view of 
1 
2 



J (I6(4)l 2 - I6(4)l 2 ) =e x -e 2 - v(^(4)) + v(z 2 (4)) 

= E 1 -E 2 - V'(x*) [^(4) - z 2 (4)] 



with a;* S [x2(t' 1 ), we have 

|VV) M4) - x 2 (4)] | < e ct [z x (4) - z 2 (4)] < e 7 e ct . 
Therefore, if e 7 e c * < 1, we have E± — E 2 > and 



1 



16(4)1- I6(^| 2 >^(Ei-E: 



1 



2 

The same holds for t' 2 , which yields 

(6(4) - 6(4)) - (6(4) - 6(4)) > ^ c \e y - e 2 ), 

whence the existence of a constant c > such that 

14 - 4| ^ ce- 2CT {E 1 - E 2 ) and \j' e \ ^ ce- 2CT (£i - E 2 ). 
The number N e (t) of intervals of the type J' £ satisfies 

N £ (t) x ce- 2Ct {E 1 -E 2 ) < t 
whence the second point of the claim. □ 
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